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'\jr{ ■ Abstract 
O 

■ We compare the lattice results on QCD phase diagram for two and three flavors with the hadron 
^ . 

. resonance gas model (HRGM) calculations. Lines of constant energy density e have been deter- 
(D ■ 

pLn ! mined at different baryo-chemical potentials fj^B- For the strangeness chemical potentials //5, we 



use two models. In one model, we explicitly set /ig = for all temperatures and baryo-chemical 



. potentials. This assignment is used in lattice calculations. In the other model, fis is calculated 

> : 

, in dependence on T and fiB according to the condition of vanishing strangeness. We also derive 

cn 

' an analytical expression for the dependence of Tc on ^s/T by applying Taylor expansion of e. In 



both cases, we compare HRGM results on Tc — fis diagram with the lattice calculations. The 
agreement is excellent, especially when the trigonometric function of e is truncated up to the same 
Q^l order as done in lattice simulations. For studying the efficiency of the truncated Taylor expansion. 



we calculate the radius of convergence. For zero- and second-order radii, the agreement with lat- 
tice is convincing. Furthermore, we make predictions for QCD phase diagram for non-truncated 
expressions and physical masses. These predictions are to be confirmed by heavy-ion experiments 



X 

' and future lattice calculations with very small lattice spacing and physical quark masses 
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I. INTRODUCTION 



The QCD phase diagram at finite temperatures and densities attracted increasing at- 
tention [l|, especially as it became possible to perform lattice QCD simulations at finite 
baryo-chemical potential hb 

. The numerical studies for the equation of state at 
finite provides a valuable framework for understanding the experimental signatures for 
the phase transition from confined hadrons to quark-gluon plasma (QGP). The heavy- ion 
experiments are aiming to explore the QCD phase diagram. Therefore, it is of great inter- 
est to show the interrelation between strangeness and baryo-chemical potentials and T^. in 
hadron resonance gas model (HRGM) compared with the available lattice QCD simulations. 

It is known that QCD phase diagram has a very rich structure. From numerical simula- 
tions, we know that the location of the phase transition line depends on the quark masses 
and fiavors and the way of including the strange quark chemical potential /i^ at different 
values oi ^B and T. The isospin chemical potential can play an additional role. But relative 
to fiB and fis, the isospin chemical potential is very small, so that we can assume an entire 
symmetry in light quark potentials and therefore ignore the isospin chemical potential. 

The first point in QCD phase diagram, namely the point at T^. and u = 0, has been a 
subject of d«.e,. >att,ce s,.u.at,„„s flflflflBBa. We .now so fa. that fo. two ,ua* 
fiavors {uf = 2) the transition is second order or rapid crossover and the critical temperature 
is Tc ~ 173 ±8 MeV. For = 3, we have a first order phase transition and ~ 154 ±8 MeV. 
For Uf = 2-1-1, i.e. two degenerate light quarks and one heavy strange quark, the transition 
is crossover and ~ 173 ± 8 MeV. For the pure gauge theory, T^. ^ 271 ± 2 MeV and the 
deconfinement phase transition is first order. 

The lattice QCD simulations at /i^ 7^ is still lacking an effective exact algorithm 
and suffer from the sign-problem. The fermion determinant gets complex and therefore the 
conventional Monte Carlo techniques are no longer applicable, since the lattice configurations 
can no longer be generated with the probability of the Boltzmann weight. However, during 
the last few years considerable progress has been made to overcome these problems 0, 0, 0, |J 

The significant numerical results on positioning the QCD phase diagram we have so far 
is that the transition line can be described by a parabola. This simple relation can be 
viewed as a refiection of the truncations done in the Taylor expansions of different thermo- 
dynamic quantities calculated on lattice. In addition, we know from effective models such 
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as bootstrap and Nambu-Jona-Lasinio models that the structure of the phase diagram is 
complex. The freeze-out curve takes a much different behavior at large chemical poten- 
tial 0, Q| . Nevertheless, one might think that for small chemical potential {^q ~ T^) the 
curvature of Tc-dependence upon fig can be fitted as a parabola, where fig is the quark 
baryo-chemical potential. The situation at very large chemical potentials is not clear. One 
might need to take into account other effects, such as quantum effects at low tempera- 



tures 
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16l |. which might be able to describe the change in the correlations 



from confined hadrons to coupled quark-pairs. 

In present work, we take advantage of our previous work [12l,llSl,ll2i on analysis the critical 
temperatures for different quark masses and on reproducing the lattice thermodynamics at 
zero and finite fiB^Y HRGM. We assumed that the deconfinement is driven by a constant 
energy density. We have shown that the degrees of freedom rapidly increase at 

eATr ^ Ur^ = 0). Here, we assume that remains constant along the whole phase transition 
line [19, 20], ec(Tc,0) = edTcUB)- Concretely, we propose that the existence of different 
transitions does not affect the assumption that ec is constant for all /i^-values. We have to 
emphasize here that it is not possible in the framework of this model to make any statement 
about the transition at very large fis and low T. As T ^ Tc, HRGM is no longer applicable. 
The nature of the degrees of freedom in this region is very different from that of the nearly 
non-interacting QGP at high T and low /i^. The message we have in this paper is that the 
condition driving the QCD phase transition at finite T and fiB ll^ is the energy density. 

Its value is not affected by the conjecture of existing of different transitions along the whole 

fiB-3X\S. 

The model will be presented in Sec. |n] In Sec. IIIH we introduce expressions for the 
lines of constant Tc. The lines of constant energy density are given in Sec. IIVI The results 
are discussed in Sec. El Sec. EH is devoted to the radius of convergence. In Sec. IVIIl we 
summarize the conclusions. 



II. THE MODEL 

Assuming an ideal quantum gas consisting of point-like hadron resonances, the canonical 
partition function for one particle and its anti-particle reads 
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V r°° 

Z{V,T,fi) = g-- rfA;A;2{ln(l±e-("-^)/^) + ln(l±e"(^+'^)/^)}, (1) 
2vr^ Jq 

where ± stand for bosons and fermions, respectively, e = (k"^ + m^)^/^ is the single-particle 
energy and g is the spin-isospin degeneracy factor. Under the given assumptions, one can 
sum up the contributions from all resonances pieces, so that 

oo 

\nZ^''\V,T,fi) = J2^nZ,{V,T,fi,). (2) 

i 

In this expression, there are two important features included; the kinetic energies and the 
summation over all degrees of freedom and energies of resonances. On the other hand, 
we know that the formation of resonances can only be materialized through strong inter- 
actions j^l; Resonances (fireballs) are composed of further resonances (fireballs), which in 
turn consist of resonances (fireballs ) and so on. 

In spite of this, if one would like to take into consideration all kinds of resonance inter- 
actions in HRGM, then, for instance, by means of the S'-matrix, we can re-write Eq. (j2)) as 
an expansion of the fugacity term. 

oo 



u=2 



S-matrix describes the scattering processes in the thermodynamical system j22]. a,y{T) are 
the so-called virial coefficients and the subscript v refers to the order of the multiple-particle 
interactions. 

n r°° f) 
aAT) = i-J d^ 5^(2/ + (4) 

The sum runs over the spatial waves. The phase shift 5i{w) of two-body inelastic interactions, 
for instance, depends on the resonance half-width F^, spin and mass of produced resonances. 



lnZ(-)(T/,T,/.) = lnZ(^^)(T/,T,;.) + |^ / dw (5) 

By inserting — /i in place of /i in Eq. (jSl), we take into consideration the two-particle inelastic 
interactions, from which the anti-particles will be produced. For narrow width and/or at low 
T, the virial term reduces, so that we will get the non-relativistic ideal partition function of 
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Fig. 1: The strange quark chemical potential /Ug vs. T/Tc for [iqjT = 1 and fiq/T = 0.5 (straight 
lines). The results are fitted according to Eq. Q. At T = 0, we find that fiq = fj,s = 0. As T — > Tc, 
the strangeness chemical potential approaches the baryo-chemical potential, //s = f^q- The units 
used here are ^/a ~ 420 MeV. 



the hadron resonances with effective masses M^,. In other words, the resonance contributions 
to the partition function are the same as that of free particles with some effective mass. 
At temperatures comparable to F^, the effective mass approaches the physical one. Thus, 
at high temperatures, the strong interactions are taken into consideration via including heavy 
resonances, Eq. (EJ. We therefore suggest to use the canonical partition function Eq. (EJ without 
any corrections. Furthermore, we do not apply the excluded volume corrections. We include 
all hadron resonances with masses up to 2 GeV, such a way we avoid the singularities 
expected at Hagedorn temperature 13| • 

In the two sections which follow, we discuss how to include the strangeness chemical 
potential f^s in HRGM and lattice QCD. 



A. Us in hadron resonance gas model 



The /iac^ron-based chemical potentials fiB, fJ'S are related to the quark-hased ones, fig, fis 



(6) 
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Assuming that the isospin and charge chemical potentials are vanishing, we use the following 
combination for the hadron resonances 



11 = Sbfig + sfis, (7) 

where b and s are the baryon and strange quantum numbers, respectively. Obviously, this 
expression is valid for baryons as well as for mesons. The quantum numbers are entirely 
conserved. 

The initial conditions in heavy-ion collisions apparently include zero net strangeness. 
This is expected to remain the case during the whole interaction unless an asymmetry in 
the production of strange particles happens during the hadronization. As we are interested 
in the hadron thermodynamics and location of QCD phase transition, we suppose that the 
net strangeness is entirely vanishing. The average strange particle number reads 



N 

i ■ 

< n, > 



1 ^ (,)91nZ«(V^ 



\nZ is given in Eq. Q and fis is given in Eq. ©• = exp(/is/T) is the fugacity factor of 
s quark. The procedure used to calculate the quark-hased fig is the following: For given T 
and fiq (or yUs), we iteratively increase /i^ and in each iteration, we calculate the difference 
< Us > — < rig >, Eq. (jHl). The value of fig which disposes zero net strangeness is the one 
we read out and shall use in calculating the thermodynamic quantities. As in Eq. © the 
relation between fig and fis is given by taking into consideration the baryonic property of 
s quark. The resulting fig for different fig (or /i^) and T are depicted in Fig. [TJ From this 
numerical method, we fit function of T and /is, 

0.138 1?^=^ 
^ 1-2.4 ^^2 + 2.7^^3' 

where = fJ^q/T and 6 = T /T^.. 

Let us note here that for these calculations we have re-scaled the resonance masses in order 
to be comparable to the quark masses used in lattice QCD simulations. The procedure 
of giving unphysically heavy masses to the hadron resonances is introduced in 0, isl ]. 
In these calculations, the mass of lightest Goldstone meson becomes = 770 MeV and 
consequently, the critical temperature gets almost as large as 200 MeV at iib = 0. 

In order to guarantee vanishing net strange particle numbers in HRGM as the case in 
heavy- ion collisions, it is not enough to simply set /Xs = and consequently, f^s = f^g = I^b/'^ 



in Eq. (jH)). However, there are publications in which the authors have assigned /is to zero in 
hadron matter and afterwards apphed the Gibbs condition for the first order phase transi- 
tion to QGP. The reason is obvious. Aside the baryons, the strange mesons with different 
contents of s quarks play determining roles at different temperatures and therefore, affect 
the final results, Eq. Setting /i^ = 0, leads to violating the strange quantum numbers. 
Nevertheless, we will show here calculations in which we set fig = 0. We do this in order 
to extensively compare with lattice results. We apply this assignment in order to check the 
ability of HRGM in reproducing the current lattice simulations [2^. After accomplishing 
this successfully, we can go beyond the lattice constrains to show the physical picture. We 
will make predictions for QCD phase diagram for physical masses and non-truncated Taylor series 
for thermodynamical expressions. These predictions are to be confirmed by heavy-ion experiments 
and future lattice calculations with very small lattice spacing and physical quark masses. 

For completeness of the discussion, we recall the situation in the plasma regime (see also 
next section). For conserving strangeness at T > Tc, we have to suppose that fis = 0. 
consists of one baryonic part /i^/S and another part coming from the strangeness quantum 
number —^s- From Eq. (jH}, we then get 



This result is numerically confirmed in Fig.^ For /ig = (or = 0), we find that /is = for 
all temperatures, /is increases with increasing both /ig and T. At Tc, we find that fiq ~ /is. 
Therefore, we can suggest to set fig = fis for all temperatures above T^,. 

We can so far summarize that /is in the hadron matter has to be calculated in dependence 
on fiB and T under the assumption that the net strangeness is vanishing. In the QGP 
phase, one might fulfill this assumption by setting /is = /is/S. As we will see later, in lattice 
simulations, one assigns /i^ = /is = 0. We deal with all these cases in this work. 

B. fig in lattice QCD 

In the Euclidian path integral formulation, the partition function of lattice QCD at finite 
T and /i reads 



f^q — /i_B/3. 



(10) 



Z(T,/i) = Tre-(^-^^)/^ 




(11) 
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where i^) and A are the fermion and gauge fields, respectively. The chemical potential fi 
is given in Eq. ((7j). By Legendre transformation of the Hamiltonian H, we get the Euclidian 
action S = J^^^ dt jy d?xC The fermionic action is 



1 4 _ 



2 

k=l 



(12) 



where a is the lattice spacing. As given in Eq. (jH}, the number density of the quarks with 
flavor number x is obtained by derivation with respect to 

d 

fix = -^\nZ{T,ij,x). (13) 

For checking the dependence of fig on /i^ and consequently on /i^, it is enough to approximate 
the fermionic part of lattice QCD Lagrangian for three quark flavors as 

\xe{u,d} j 
^ figUu + figUd + HsUs- (14) 

To taken into account the conservation of the baryon and strange quantum numbers, the 
summation in last expression has to run over s quarks, too. By doing this, last term turns 
to be {fiq — fig) Us. Then we expect that the strangeness on lattice vanishes at fig = fJ'q- But 
from Eq. ()14|) . which reflects the situation in current lattice simulations, we flnd that = 
for fig = 0. 

As discussed in previous section, the strangeness in QGP is conserved ai fig = fig = fis/S. 
In the hadron regime, especially at large /i^, /is (or fis) has to be calculated as a function of 
T and fiB (Fig-C])- In spite of these considerations, the reliable lattice QCD simulations are 
still limited to ~ 3Tc. As we will see in Sec. |3 at this small value, there is practically 
no big difference between fig = and fig = f(T,fiB)- 
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III. LINES OF CONSTANT 



A. Tc{p) in hadron resonance gas model 



In Boltzmann limit, the energy density in an ideal quantum system consisting of one 
particle and its anti-particle can be expressed as 



cosh ( 

T 



(15) 



fi is given in Eg. (I7|l . We can divide this expression into two sectors: one meson- m and one 
baryon-sector b [3 Q , 



ef,{T,fi) = e^{T)+e,{T) cosh(/i/T). 



(16) 



Applying truncations in the Taylor expansions of cosh (/i/T) up to the second order of 
/i/T - as done in lattice calculations j2J] - we can estimate the lines of constant in 
HRGM. In doing this, we assume that eh(T,fi) = eh(Tc,0). Furthermore, we assume that 
the dependence of energy density on the quark mass is quite small. We see this in Fig. |21 
which is reported in Ref. [2^. At Tc and for the same temporal lattice dimension A^,-, there 
is very small change in Sc for different quark masses. The latter are related to the pion 
masses via oc rriq. 

The starting point in expressing Tc^Hb) for constant ratios UbIT is to expand Eq. (fTBj) 
around the points T = Tc and fiB = 0- The second-order expansion gives 



eh{T,fXB) = £h(T„0) + 
= eh{Tc,0) + 



dT ^ ^' 2^^ 



dT 



(T - Tj + i5,(T„ 0) [^y 



Under the assumptions given above, this leads to the following parabola: 



B) 



Te(/iB = 0) 



1-^ 



2 Tcii^B = 0) 



gb(Te, 0) 



dT 



tJ 



(17) 



(18) 



where fig = /i^/S. In order to map out the QCD phase diagram by using this analytic 
expression, we merely need to calculate the baryon energy density e^, at and fig = and 
the derivative of the hadron energy density with respect to T, deh(Tc, 0)/dT. The results are 
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given in Sec. El Evidently, it is possible to extend the above expression to include further 
higher terms of ^iq/T. 

Then to determine at given /i^, we apply besides the above analytical method the 
condition of constant critical energy density jiol. I20I (see Sec. I1V|) . In this case, T^. is defined 
as the temperature at which the energy density in HRGM reaches a certain critical value. This 
value can be taken from lattice QCD simulations at hb = 0. 



B. Tciij) in lattice QCD 

In lattice QCD simulations, the critical temperature Tc(/i) is to be calculated from the 
pseudo-critical coupling /?c(At) by determining the susceptibility peak in either the Polyakov 
loop or the chiral condensate in (3 — ^ dimensions. The lattice beta function (3{a) which 
can be obtained from the string tension is needed in order to express the results in physical 
units. From the first non-trivial Taylor coefficients of Tc(yu), we get 

where N^- is the temporal lattice dimension. 

n 

In the lattice QCD simulations |^ with Uf = 2 and quark mass aruq = 0.1, the dependence 
of Tc on fig has been found to take the following (parabola) expression: 

'^"^^'^ - 1 - 0.070(35) ( ' . (20) 



Tc{l^q = 0) ' ^V^c(/ig = 0), 

In other lattice QCD simulations with the same flavor number but quark masses four times 
heavier than the physical masses j^. 



2 



1 - 0.050(34) . (21) 

Tc(/x, = 0) ^ ^ V^c(/i, = 0)y' ^ > 

It has been concluded that the last relation remains almost unchanged for aniq = 0.005. 
In this regard, we have to remember that for small quark masses perturbative beta function 
has been used. As we will see later, we actually flnd an increase in the curvature with 
reducing resonance masses from the values which very well simulate the current lattice 
calculations (m,r ~ 770 MeV) to the physical masses, at which the lightest Goldstone meson 
is the physical pion (m^ = 140 MeV). 
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The results from Eq. (pUjl are depicted in Fig. El In the same figure, we plot Eq. (j^D). We 
also plot other lattice results j^, namely the short vertical lines which give Tc at constant 
energy density. 

I I 

Uf = 3 lattice QCD results ^6] for different quark masses obtained so far can be sum- 
marized as 



- 0.025,6, (^)^, .»,^0.1, (22) 

= 1 - 0.114(46) ( y^^^'^Q) ) , anig = 0.005. (23) 

Again, in these simulations the beta function for the small quark masses has been calculated 
perturbatively. The results from these two expressions are depicted in Fig. El and Fig. 
respectively. From other n/ = 3 lattice QCD simulations j^], 

J^J^^ = 1 - 0.0610(90) (I)' + 0.00235(89) (^)\ (24) 

The curvature from Uf = 2 + 1 lattice simulations [3| has been compared with the three 
fiavor one ^ and concluded that there is a complete agreement. As we will see in Fig. El the 
most recent Uf = 2 + 1 lattice simulations of j2^ result in a curvature smaller than those 
from other lattice simulations. This can be understood according to the different fermionic 
actions used. 



IV. LINES OF CONSTANT ENERGY DENSITY 

In previous work we have used HRGM in order to determine Tc corresponding to 
a wide range of quark (pion) masses at fiB = 0. The masses range from the chiral to pure 
gauge limits. We have seen that the condition of constant energy density can excellently 
reproduce the critical temperature Tc as a function of and nj. In present work, we 
extend this condition to finite chemical potentials fiB- We use two models for including the 
strangeness chemical potential /ig. In the first model, we calculate /i^ in dependence on /i^ 
and T according to the condition of zero net strangeness. In the other one, we explicitly 
set the quark-based fis = 0. The last case is used in lattice QCD simulations. As given in 
Sec. IIIBl the proper inclusion of fig in QGP and in lattice Lagrangian is to assign to it the 
same value of iib/^- Nevertheless, both choices can be accepted for the current lattice QCD 
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Fig. 2: The energy density normalized to given in dependence on T/Tc for different quark masses 
and temporal lattice dimensions Nr- This figure is reported in Ref. j25l] . 

simulations, since the most reliable lattice calculations are currently performed at small 
baryo-cliemical potentials (/ig ~ Tc). Consequently, the strangeness chemical potential is 
expected to be much smaller or at lest as small as the baryo-chemical potential. 

As discussed above, in HRGM, the energy density at finite chemical potential can be 
divided into two parts: one from the meson sector and another one from the baryon sector. 
For the first part, we can completely drop out the fugacity term. For symmetric numbers of 
light quarks, the baryo-chemical potential of mesons is vanishing. But for strange mesons 
the strangeness chemical potential assigned to their s quarks should be taken into account. 
For the baryon sector, the chemical potential is given by Eq. ((Tj). 

The question we intent to answer is: which value has to be assigned to the critical energy 
density at finite chemical potentials? We recall the lattice QCD simulations. In Fig. |2l we 
see that the energy density j^^] at Tc is not a singular function but can rather be defined 
at different critical values depending on and nig. This reflects the nature of the phase 
transition, cross-over. On the other hand, the uncertainty in calculating the energy density 
on lattice is very large, almost a factor 2. The reason for this is the uncertainty in estimating 
Tc, Eq. (pUj) . As mentioned above, Tc is determined according to maximum susceptibility. 
The uncertainty in this quasi Tc is ~ 10%. The coefficients of T have additional ~ 20%. 

The critical energy density at which we define Tc{^b) in HRGM is taken from lattice 
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QCD simulations at /i^ = '8|. In lattice units, the dimensionless energy densities for 
Uf = 2 and 2 + 1 are e/T^It^ = 4.5 ± 2 and = 6.5 ± 2, respectively. We take an average 
value and express it in physical units. Thus, ec = 600 ± 300 MeV/fm^. 

As in [4], we assume that this value remains constant along the phase transition line; fis- 
axis. The existence of different phase transitions (cross-over and first-order) and the critical 
endpoint, at which the transition is second-order, is assumed not to affect this assumption. 
As mentioned in Sect. IIII Bl the phase transition line on lattice is to be estimated according 
Eq. (fT^ . i.e. up to (n/T)^. Additional to this method, the hne of constant energy density 
has been calculated [J]. By comparing the two lines, it has been found that they are almost 
coincident. In other words, the constant ec can be used to estimate the critical line at finite 
fiB- We have to remember that the lattice estimation of e has a very large uncertainty, 
e oc T^. Therefore, the transition line according to constant ec has a larger uncertainty, 
because of the additional uncertainty in the derivatives of ec with respect to hb and /x^. 

What are the consequences if the assumption of constant ec(/iB) disregarding the un- 
certainty turned out to be incorrect? To discuss such consequences, we first recall the 
second-law of thermodynamics. 

de = T ds — p + hb duB, (25) 

where s, p and hb are the entropy, the pressure and the baryon number density, respectively. 

If ec were a decreasing function of /x^, this means that low incident energies are much more 
suitable to produce the phase transition from confined hadrons to deconfined QGP than the 



3, i.e., 



high incident energies! But the low limit should be given by the freeze-out curve 
the hadronization phase diagram. We know from phenomenological observations that both 
freeze-out curve and phase transition line are coincident at low /xr (high incident energy). 
For large fiB the two lines are separated. The freeze-out curve 10] is given by s/T^ = 7, 
entropy driven. Then along freeze-out curve it is expected that e slightly increase with fiB- 
This means the assumption that ec decreases with increasing fiB will give a phase transition 
which has much smaller than the freeze-out temperature. 

In the other case, that ec increases with increasing /ig, we expect the ec required for the 
phase transition gets larger with decreasing the incident energy! According to Eq. (jSHI), 
the phase diagram is given by = de/ds at constant ub and ^b = de/dnB at constant 
T. Then for increasing ec{fiB), the critical temperature is expected to increase or at least 
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remain constant. This result has the conconsequences that the phase transition at large 
fiB will be diffecult to materialize in heavy- ion collisions or impossible. Also the phases 
of coupled quark-pairs (color superconductivity) will be expected for very large /i_B or not 
allowed at all. 

Since we use constant Cc down to low temperatures, one might ask if there could be any 
relation between ec and the experimental value of nuclear density. Obvioulsy at T = 0, 
HRGM, Eq. is no longer applicable. In this limit, we suppose that the hadron gas is 
composed of degenerate Fermi gas of nucleons. We can therefore calculate fiB corresponding 
to the normal nuclear density at T = 1^. The value is 979 MeV. The energy density in 
this limit reads 



e(/i 



B 



4 

m 



4-^-2 



m 



- 1 



/4 



In 



+ 



(26) 



V. THE RESULTS 

Since we want to compare HRGM with lattice results 0, 0], some details about the 
lattice simulations at fig ^ are in order. Extensive details about HRGM are presented 



in Refs. 



HQ, 



19|. In order to avoid the sign problem in lattice simulations [4|, the 



derivatives of thermodynamic quantities with respect to /i^ = at the point Tc are first 
computed and then their Taylor expansion coefficients in terms of finite fig are calculated. 
The curvatures according to Tc{^b) and ec{fiB) are derived for different rriq. For instance 
for arriq = 0.1, the corresponding pion mass in lattice units for = 2 is am.^ = 0.958(2), 
where aa = 0.271(10) Q]. This leads to = 773 MeV. Keeping these features in mind, 
we performed our calculations for physical and re-scaled resonance masses. In Fig. El we 
find that our results are coincide with the lattice ones. In Ref. [3], the lattice simulations 
are performed with an imaginary chemical potential for Uf = 2 staggered quarks. For 
fi = ifij, the sign problem is obviously no longer existing. Using analytic continuation of 
the truncated Taylor series, one can then go to real fi. The light quark mass is four times 
the physical mass. The location of at /i^ = has been extrapolated to the chiral limit. 
The lattice results are represented by the short curves in Figs. El IH and IHl 
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A. Results for two flavors 



From the canonical partition function Eq. (^J), we can derive the energy density at finite 
chemical potential yU 7^ as 

, BT\nZ{T,^i) , dT\nZ(T,fi) 
e(T,/x) = T -±^-T\nZ{T,fi)+f, 

^ k'^ dk _L 1 • (27) 



27r2 Jo el^e^)-/^]/^ ± 1 ■ 

In Boltzmann limit and by taking into consideration only one particle and its anti-particle, 
we get the expression given in Eq. (fT3j) . It is obvious that the trigonometric function included 
in last expression are not truncated. In calculating this quantity in HRGM, we sum up over 
all resonances we take into account. 

We first start with two flavors. It is not needed to care about the strangeness chemical 
potential. As discussed above, there are two lattice QCD simulations for nj = 2. In the 
first one, re-weighting methods, the physical quantities which usually can be calculated at 
^ wthout a,>y difflculfes have responses at fin.te The responses are ntiUzed to 
estimate the thermodynamic quantities at finite fiB- The other lattice simulations p| use 
imaginary fiB and afterwards apply analytical continuation. At given fis, Tc is determined 
as the temperature at which the energy density e equals 600 MeV/fm^. 

The results are plotted in Fig. El The solid circles represent the results for re-scaled 
resonance masses. The open circles represent the results for physical masses. The two lines 
connecting the points are obtained by fits according to: 
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Te(/.. = 0) ' '^U(/^, = 0)J • ^'^^ 

The /iq- values are restricted within the range < /ig < Tc. /ig is given in Eq. For the 
re-scaled heavy masses, the fit parameters are Ci = 0.115(36) and Tc{fiB = 0) = 196.3 MeV. 
Plugging these parameters in Eq. (j^H|) . we get the top line in Fig. El The comparison 
with the lattice calculations Eq. ()20|) gives a satisfactory agreement, especially at low /i^. 
Nevertheless, it is obvious that our results at large fiB he below the lattice ones. The reason 
for this discrepancy will be discussed later. 

One might ask whether the pion gas with very heavy masses would be able to describe the 
lattice results? To answer this question, we refer to ITj, ll8[ . In order to simulate the lattice 
QCD thermodynamiscs, for instance the rapid increase in e near T^, we need to include the 
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Fig. 3: T — fiB phase diagram for two quark flavors nj = 2. The vertical lines give the lattice results ^. 
The short lines show the results according to a constant ec- The long ones are for constant Tc- The Lattice 
simulations are performed for large quark mass. The corresponding Goldstone pion gets a mass of 770 MeV. 
The solid circles give our results heavy quark masses. The results for physical quark masses are given by 
the open circles. 

corresponding degrees of freedom in HRGM. A hadron gas of pions can reproduce ~ 15% 
of the lattice QCD thermodynamiscs. Including heavier resonances gives the correct QCD 
thermodynamics. Furthermore, for lattice calculations with very heavy quark masses (pure 
gauge) we need to include the low-lying glueballs |l^. 




As given above, the coefficient in the front of {fiq/TY in Eq. ()18|) can be calculated in 
HRGM. For non-strange hadron resonances with re-scaled masses, the coefficient gets the 
value 0.0767. Comparing to Eq. (j^UI) . this value is obviously much better than the value of 
ci in describing the lattice results. As we will see, this discrepancy is to be related to the fact 
that ci has been calculated from fitting Tc at finite /Ub in the results which we have obtained 
from non-truncated e{T,fiB)- We should emphasize that all Taylor terms that are higher 
than the second one are explicitly excluded in deriving Eq. ()18|) . This partially explains the 
disagreement between HRGM and lattice results shown so far. We also study the case for 
physical resonance masses. The fit parameters are Ci = 0.195(21) and Tc{fi = 0) = 175 MeV. 
The coefficient in Eq. p8|l is found to be 0.1368. 

In the following, we confront the lattice results with HRGM results which have been 
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Fig. 4: The same as in Fig.O Here, we show the results obtained from the condition of constant truncated 
up the different Taylor orders. The truncation is obviously able to describe the lattice QCD simulations for 
different quark masses ^ 22|. The solid lines give our predictions for non-truncated trigonometric functions. 



obtained from truncated expressions. We show the results in Fig. |3J We start with results 
from entire Taylor expansion for heavy resonance masses (solid line). This is equivalent to 
use Eq. ()15|) instead of Eq. ()27|) . We note that this line shows almost the same behavior 
as that of the solid circles in Fig. IHl In deriving Eq. (fTH|l . the Boltzmann limit is assumed. 
The results from different truncations are also plotted. We find that the energy density 
truncated up to the second or fourth order of /i/T can describe the lattice results better 
than the solid line. The agreement between the curvatures of these two lines and that from 
the analytical expression Eq. (fTHj) is excellent. The ability of truncated expression to produce 
results comparable to the lattice ones is to be explained by the fact that the lattice results 
themselves have been obtained from truncated thermodynamic expressions. 

We plot in the same figure the HRGM results with physical masses. The curves obtained 
from different truncation terms in Eq. (1T5|) represent our predictions when it will be possible to 
perform lattice simulations for two quark flavors with physical masses. These predictions has 
to be checked by future lattice simulations. We note that the solid line is comparable to the 
bottom points plotted in Fig. 01 The quark masses used in Ref. [^^I are relative heavy. 
Nevertheless, we see that the second order (dashed line) agrees very well with these lattice 
results j^^l (solid circles). 
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Fig. 5: The T~ hb phase diagram for riy = 3. The vertical short lines represent the lattice results 4]. The 
above band corresponds to heavy quark masses. The bottom one is for physical masses. The solid circles 
are our results for heavy resonances and for /is = 0. The open circles are for /is = f{fj,B,T). The squares 
give the results for the physical masses. The curves are fitted according to Eq. (|28|l . 

B. Results for three flavors 

We have to include strangeness chemical potential fis (Sec. Ill A|l . Assuming that the 
three quarks are degenerate, we can use the same critical energy density as done in previous 
section. The results are shown in Fig. El The solid circles show the results for re-scaled 
heavy resonance masses. As done in lattice calculations, we first set fis = 0. The resulting 
points (solid circles) are fitted within the range < fig < according to Eq. (j^H|) . The fit 
parameters are Tc = 190 MeV and Ci = 0.1016. The coefficient in front of (/x/T)^ in Eq. (fTH)) 
is 0.0505. Comparing with Eq. ()22|1 . this value can describe the lattice curvature much better 
than ci. We also calculate /i^ in dependence on fiB and T assuming strangeness conservation. 
The open circles show the T^ — diagram corresponding to this value of yU^- We note that 
the former case (/x^ = 0) is much closer to the lattice results than the latter one. The reason 
is, as we mentioned above, that in lattice calculations used to be assigned to zero. 

The results for physical masses are also shown in Fig. El The solid squares represent 
the results at /i^ = 0. These points are fitted according Eq. (j2Hl)- The fit parameters are 
Tc = 164 MeV and ci = 0.17. The coefficient in Eq. ^ takes the value 0.1122, which 
agrees very well with Eq. (j^ . The open squares represent the results at jis that depends on 
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Fig. 6: T — phase diagram as in Fig.|Sl Here, we check the effects of truncated trigonometric functions. 
The top panel shows the results for heavy masses. Truncating e(r, /is) up to second order reproduces the 



lattice results 4]. In the bottom panel we show the results with physical masses, the lattice results 
are very well reproduced by the condition of constant truncated e. The results from non-truncated e at 
Ha = and Us = f(T, iib) are also plotted, the bottom curves, respectively. The agreement with the lattice 



simulations l27 



29| is also convincing 



fiB and T. The results with re-scaled resonance masses are given as circles. For fi^ = 0, we 
get results closer to lattice results than for jig = f{T, hb)- In this figure, the energy density 
is calculated according to Eq. (j^Tj). In other words, the results plotted here are deduced 
from expression like Eq. ()15p. i.e. without any truncations. 

We plot in Fig. EJ the results from different truncations, in order to compare HRGM 
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results with lattice simulations. For heavy masses the results are given in the top panel. 
The results for physical masses are given in the bottom panel. It is clear that the truncation 
up to the second order gives results in a good agreement with the lattice simulations |^|. With 
higher truncations, we get curves with larger curvatures. We also plot two curves from 
non-truncated expressions for the trigonometric functions. The solid curve is obtained by 
computing fis in dependence on T and /i^- Obviously, it is identical to the solid symbols 
plotted in Fig. The dashed curve - next to the solid one - shows the results in which /i^ 
is entirely vanishing. It gives the same behavior as that of the open symbols in Fig. El 

We can so far conclude that the lattice results [4J can be reproduced by HRGM, if the 
trigonometric functions are truncated in the same way. The lattice curvature within the region 
< < Tc can excellently be described by HRGM. 

In the bottom panel, we show other three flavor lattice results. In Ref. j|27|, the numerical 
simulations have been performed with Wilson gauge action and three degenerate flavors of 
staggered fermions. The quark masses are ranging between 0.025 < aruq < 0.04. It is clear 
that our results can reproduce the structure given by these data. We see that the curva- 
ture can be described by second or fourth order in HRGM. This is also valid for n f = 2 + 1 

n 

lattice results [29||. In this case, the quark masses arriu^d = 0.0092 and am^ = 0.25. Chiral 
extrapolations are done by heavier light quark masses. We also plot the latest calcula- 
tions of the location of the endpoint. The endpoint in [s^ lies at /i^ = 420 MeV. The 
corresponding temperature has not yet been calculated. The endpoint in [2^ has the co- 
ordinates /i_B = 360 ± 40 MeV and T = 162 ± 2 MeV. As discussed above, we assume that 
the existence of endpoint does not affect our results. 



VI. RADIUS OF CONVERGENCE 



As we have seen, for a reliable comparison with the current lattice simulations, we have 
to apply truncations in the thermodynamic expressions in HRGM. Our objective here is 
to check the efficiency of truncated series in locating the phase diagram. The radius of 
convergence reflects the singularity near T^. It approaches unity near T^. The energy den- 
sity normalized to can be expressed as a trigonometric function depending on the ratio 
fi/T fSec.lIIl|). We use the property that the energy density is an even function in fJ^s/T 
and therefore write 
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Fig. 7: Comparison between radii of convergence from the lattice |30| (points) and HRGM (lines). There 
is an excellent agreement with the existing lattice results. The Tc-curvature has been calculated by using 
lattice results with two flavors (Fig.lSJ. The short vertical lines give the values of corresponding radii at Tc. 
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where = (Tc/T)"' 3"'/n\ and n is an even positive integer. The radius of convergence 
of Taylor expansion of the partition function, Eq. is to be estimated by the ratios of 
subsequent expantion coefficients, 
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(30) 



Since the expansion is an even series of fis/T, the square root is expected to arise. In order 
to calculate the zero order radius, we recall the results at /i = (not included in Eq. (jSH)))- 



. ^0 
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C2 



1/2 



T 



+ 1 



1/2 



(31) 



[9 Vefe(T,/i, = 0) 

It is clear that po? in constast to the other radii, is T-dependent In Fig. [7[ we find that 
close to Tc, Po = 0.982, p2 = 1.18, p4 = 1.82 and pe = 2.5. These values are shown as 
short vertical lines. The results agree very well with the available lattice results j3{|; as 
Po ^ 1. The radius of convergence approximately gives the lower bound for the critical 
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endpoint {fig > Tc). Obviously this value is higher than the recent lattice results 
In principle, one expects that due to the absence of critical behavior in HRGM, the Taylor 
expansion, Eq. (|29|) . has an infinite convergence radius for all temperatures. At and 
for the convergence radius p = 1, one expects that the hadron degrees of freedom are 
indistinguishable from the degrees of freedom in QGP. On the other hand, the radii p from 
lattice simulations are bounded from above by the location of phase transition. Therefore, 
it is expected that the radii in Eq. ()30|) stay close to unity near Tc. In fact, this is the 
case for our low-order coefficients. The agreement between our results and the lattice ones 
is convincing. There are no published lattice results on the G^'^-order radius of divergence. 
Nevertheless, according to our expansion coefficients, we expect that p are steadily changing. 

VII. CONCLUSIONS 
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Fig. 8: Summary of our results on the QCD phase transition. In calculating both curves, we use all 
hadron resonances and suppose that the isospin chemical potential is vanishing. To obtain the top curve, 
we assigned fi^ to zero, fis is calculated in dependence on T and /i^. The results are given by the bottom 
curve. We note that both curves will cross the abscissa at the point of normal nuclear density. 

We used HRGM to draw up — Pb diagram. The transition temperature from 
hadronic matter to QGP has been determined according to a condition of constant energy 
density. Its value is taken from lattice QCD simulations at zero chemical potential and 
assumed to remain constant along the entire /x^-axis. We checked the influence of s quark 
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Hs=f(T,HB) 




chemical potential /i^ on Tc- For including we applied two models. In the first one, 
we explicitly calculated /i^ in dependence on T and iib under the condition that the net 
strangeness vanishes. In the second one, we assigned, as the case in lattice QCD simulations, 
zero to /is for all T and ^b- The first case, /i^ = f{T,fiB), is of great interest for heavy- 
ion collisions. Furthermore, under this consideration, we expect that the strange quantum 
number is entirely conserved. This is not the case in the second model. Nevertheless, with the 
last assignment, the current lattice results are very well reproduced. On the other hand, one 
can apply the second model for current and future heavy-ion collisions. At BNL-RHIC and 
CERN-LHC energies, for instance, fiB (and consequently fig) is very small. We have shown 
that the proper condition that guarantees vanishing strangeness in QGP is to set fis = /^g- 
We did not check this explicitly. But it is obvious that Tc{fiB,fJ's = 0) quantitatively is not 
very much different from Tc^^ib, l^siT, ^b)) at small /i^. 

We have taken into consideration all hadron resonances given in particle data booklet 
with masses up to 2 GeV. For a reliable comparison with lattice, we have re-scaled the 
resonance masses to be comparable to the quark masses used in lattice simulations. With 
excluding the strange resonances, we compared our results with Uf = 2 lattice results. With 
including all resonances, we reproduced Uf = 2 + 1 lattice results. We note that increasing 
fiB and fig leads to monotonic decrease in Tc. The results from HRGM match very well 
with the lattice simulations, especially within the /le-range in which the lattice calculations 
are most reliable (/is/T) ^ 3Tc. The agreement turns to be excellent, when we take into 
consideration the truncations done in calculating the thermodynamical quantity e. 

Besides this excellent agreement in the indirect calculation of Tc via constant energy 
density, we found that the analytical expression of Tc(/i) up to the second order of (/u/T) 
greatly reproduced the curvatures calculated on lattice for different quark masses and flavor 
numbers. 

Fig. IHl summarizes our conclusions. The QCD phase diagram is plotted for a system 
including light as well as strange quarks. The Taylor expansion of energy density is not 
truncated. The two curves represents our predictions, when it will be possible to perform 
lattice simulations for physical masses and without the need to truncate the Taylor expan- 
sion. As T — s> 0, we note that both curves will cross the abscissa at almost one point. It is 
obvious that this point is corresponding to the normal nuclear density. The latter is related 
to fiB ~ 0.979 GeV. The nature of the phase transition at very low temperatures does not 
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lie within the scope of this work. However, there are many indications that the transition 
at T = occurs according to modification in the particle correlations. Changing the corre- 
lation leads to quantum phenomena, like quantum entropy 12, 13, 14, 15, Q|. We also 
note that by switching on fig an increase in Tc is expected. At small fiB, the two curves are 
coincide. 



Acknowledgments 

We gratefully acknowledge the useful discussions with David Blaschke, Frithjof Karsch, 
Berndt Miiller, Krzsyztof Redlich, York Schroder and Boris Tomasik. 



[1] K. Rajagopal and F. Wilczek (2000), hep-ph/0011333. 

[2] Z. Fodor and S. D. Katz, Phys. Lett. B534, 87 (2002), hep-lat/0104001. 

[3] P. de Forcrand and O. Philipsen, Nucl. Phys. B642, 290 (2002), hep-lat/0205016. 

[4] C. R. Allton et al., Phys. Rev. D66, 074507 (2002), hep-lat/0204010. 

[5] M. D'Elia and M.-P. Lombardo, Phys. Rev. D67, 014505 (2003), hep-lat/0209146. 

[6] R. V. Gavai and S. Gupta, Phys. Rev. D68, 034506 (2003), hep-lat/0303013. 

[7] F. Karsch, Lect. Notes Phys. 583, 209 (2002), hep-lat/0106019. 

[8] F. Karsch, E. Laermann, and A. Peikert, Nucl. Phys. B605, 579 (2001), hep-lat/0012023. 

[9] A. Tawfik (2004), hep-ph/0410329. 
[10] A. Tawfik (2004), hep-ph/0410392. 
[11] D. E. Miller and A. Tawfik (2003), hep-ph/0308192. 
[12] D. E. Miller and A. Tawfik (2003), hep-ph/0309139. 
[13] D. E. Miller and A. Tawfik (2003), hep-ph/0312368. 
[14] D. E. Miller and A. Tawfik, J. Phys. G30, 731 (2004), hep-ph/0402296. 
[15] S. Hamieh and A. Tawfik (2004), hep-ph/0404246. 

[16] D. E. Miller and A. Tawfik, Acta Phys. Polon. B35, 2165 (2004), hep-ph/0405175. 
[17] F. Karsch, K. Redlich, and A. Tawfik, Phys. Lett. B571, 67 (2003), hep-ph/0306208. 
[18] F. Karsch, K. Redlich, and A. Tawfik, Eur. Phys. J. C29, 549 (2003), hep-ph/0303108. 
[19] K. Redlich, F. Karsch, and A. Tawfik, J. Phys. G30, S1271 (2004), nucl-th/0404009. 



24 



[20] D. Toublan and J. B. Kogut (2004), hep-ph/0409310. 
[21] R. Hagedorn, Nuovo Cim. Suppl. 3, 147 (1965). 

[22] R. Dashen, S.-k. Ma, and H. J. Bernstein, Phys. Rev. 187, 345 (1969). 

[23] F. Karsch, Prog. Theor. Phys. Suppl. 153, 106 (2004), hcp-lat/0401031. 

[24] S. Ejiri et al., Prog. Theor. Phys. Suppl. 153, 118 (2004), hep-lat/0312006. 

[25] A. Ah Khan et al. (CP-PACS), Phys. Rev. D64, 074510 (2001), hep-lat/0103028. 

[26] F. Karsch et al., Nucl. Phys. Proc. Suppl. 129, 614 (2004), hep-lat/0309116. 

[27] P. de Forcrand and O. Philipsen, Nucl. Phys. B673, 170 (2003), hep-lat/0307020. 

[28] Z. Fodor, Nucl. Phys. A715, 319 (2003), hcp-lat/0209101. 

[29] Z. Fodor and S. D. Katz, JHEP 04, 050 (2004), hep-lat/0402006. 

[30] C. R. Allton et al., Phys. Rev. D68, 014507 (2003), hep-lat/0305007. 

[31] Comparing full QCD with pure gauge results, we get a feeling about the dependence of critical 
values on quark mass. ec/T^ seem to be different. On the other hand, taking into consideration 
the different critical temperatures, we find that the Ec are comparable with each other. 



25 



